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We reexamine Wald's no-hair theorem for global anisotropy in the brane world scenarios. We 
derive a set of sufficient conditions which must be satisfied by the brane matter and bulk metric 
so that a homogeneous and anisotropic brane asymptotically evolves to a de Sitter spacetime in 
the presence of a positive cosmological constant on the brane. We discuss the violations of these 
sufficient conditions and we show that a negative nonlocal energy density or the presence of strong 
anisotropic stress (i.e., a magnetic field) may lead the brane to collapse. We discuss the generality 
of these conditions. 
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It has been proposed that we may live confined to a 
four dimensional brane in a higher dimensional spacetime 
jlj; in this scenario the fundamental high dimensional 
Planck mass could be of the same magnitude as the elec- 
troweak scale, resolving one of the hierarchy problems in 
the current standard model of high-energy physics. 

This brane world scenario generally predicts deviations 
from standard model processes at high energies, a direct 
consequence of the possibility that gravitons propagate 
in the bulk Q . Nevertheless it remains viable given that 
there are no constraints on the gravitational force at very 
high energies (or very small distances); indeed current 
constraints are only beginning to probe sub-mm scales. 
One should also expect deviations from four dimensional 
Einstein gravity at high densities (or correspondingly 
strong gravitational fields). A first manifestation of this 
was seen when deriving the evolution equations for the 
cosmological scale factor, a, in the early universe: it was 
shown that a is driven by an additional term which is 
proportional to p 2 , where p is the energy density of the 
cosmological fluid on the brane [|| . 

The two standard test situations for studying grav- 
ity in the strongly nonlinear regime are the formation 
and stability of black holes and the evolution of glob- 
ally anisotropic models in the early universe. There 
have been attempts at constructing black hole solutions 
on the brane by extending the well-known four dimen- 
sional solutions to higher dimensions; however, they have 
been shown to be unstable 0. Progress is being made 
in numerically evolving the modified Einstein's equa- 
tions to produce a stable analog of the four dimensional 
Schwarzchild solution || . Some work has also been done 
on the effect of global shear on the brane, in particular 
on its effect on the dynamics of the scale factor during 
inflation §,§ . 

One of the only unambiguous predictions of the infla- 
tionary cosmology was put forward by Wald in 1983 ||: 
it was shown that initially expanding homogeneous cos- 
mological models isotropize in the presence of a positive 
cosmological constant. This leads to a firm prediction: if 
we detect any form of global anisotropy, then the universe 
cannot have undergone a de Sitter stage. As yet there is 



no detection of global anisotropy, and an upper bound on 
global shear was derived from the COBE four year data 
P]. In this paper we wish to reexamine Wald's result 
in the brane world scenario. The conditions we derive 
are more complicated than for the four dimensional case 
considered by Wald, but we can still assess their general- 
ity within the class of cosmological fluids currently being 
considered in cosmology. 

We will attack this problem locally, i.e., in terms of 
dynamical quantities on the brane; in doing this we are 
ignoring possible constraints which we may eventually set 
on the bulk due to our assumption on the brane (we dis- 
cuss the consequences of this approach at the end of this 
paper). Our starting point are the modified Einstein's 
equations derived by Shiromizu, Maeda, and Sasaki |lC| ] 
under the assumption of Z2 symmetry in the extra di- 
mension (for a five dimensional bulk). Following the no- 
tation of Maartens [Ol , these can be written as 
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where G 



and T^y are the four dimensional Ein- 



stein, metric, and energy-momentum tensors, respec- 
tively; A is the effective cosmological constant on the 
brane, and S^ v is a term quadratic in the energy- 
momentum tensor, defined by 
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The quantities k 2 /8tt and k 2 /8tt are the effective Newton 
constant on the brane and the fundamental Newton con- 
stant in the bulk, respectively. The quadratic term S pu 
represents the matter corrections of the effective four di- 
mensional gravity, and is significantly important at high 
energy, i.e., when \T^ V \ > A = 6k 2 /k 4 , where A is the 
brane tension. Finally, £^ v , which represents the pro- 
jection of the five dimensional bulk Weyl tensor on the 
brane, is symmetric and traceless, and transmits nonlo- 
cal gravitational degrees of freedom from the bulk to the 
brane. We do not make any speculation about the origin 
of the cosmological constant A which we assume to be 
positive. 
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Let us now consider two components of Eq. ([[]), the 
"initial-value" constraint equation, 



G^u" = A + ^T^u^u" 



(3) 



where u M is the unit normal to the spatial homogeneous 
hypersurfaces, and the Raychauduri equation, 



(4) 



Both GfmU^u" and R^u^u" can be expressed in terms 
of the three-geometry of the homogeneous hypersurfaces 
and the extrinsic curvature 0^,, = V„u M of these surfaces 
using the Gauss-Codazzi equations. For convenience we 
decompose M!/ into its trace O and trace-free part (7 M „, 



(5) 



where /i M „ = g^,, + u^u u projects orthogonal to u^, and 
we rewrite Eqs. (||) and (|j) as 



6 2 - 3A - -o^tr"" - I (3) i? + :iK-7 ; ,,,/*' 



(0) 



and 



6 - A + ie 2 = -o-^ff^ - (t^ - ~ 5jUI ,r) 



(7) 



where the dot denotes u^V^. Eq. ([?]) is thus the evo- 
lution equation for the expansion rate and er M „ is the 
shear of the timelike geodesic congruence orthogonal to 
the homogeneous hypersurfaces. The scalar curvature 
( 3 'R is given in terms of the structure-constant tensor of 
the Lie algebra of the spatial symmetry group (see [ p~2|Jl3| 
for Bianchi models). 

We now quickly review Wald's argument which pro- 
ceeds using the well known energy conditions on the 
energy-momentum tensor T^ v . These state that for any 
non-spacelike vector the following inequalities hold: 



T^f>0 and I - -g MV T j ft" > 0, (8) 

known as the weak and strong energy conditions, respec- 
tively. 

On the other hand, in standard gravity, the local 
quadratic term and the nonlocal term in Eqs. (^J) and (]?]) 
are absent, and it is also well known that for all Bianchi 
models, except Bianchi type-IX, < 0. Therefore, 

Eqs. (U) and (|7|) together with the energy conditions (^), 
yield the inequalities 



e < a - -e 2 < 0. 



(9) 



The second of these two inequalities tells us that O never 
passes through zero and if the universe is expanding at 
some time it will do it forever. Indeed, we have that 
O > (3A) 1 / 2 at all time. Moreover, the first inequality 
can be integrated and yields 



e < 



(3A) 



1/2 



tanh(t A 



(10) 



Thus is "squeezed" between the lower limit (3A) 1 / 2 
and the upper limit in Eq. (|Io|), which exponentially ap- 
proaches (3A) 1 / 2 on a time scale (3/A) 1 / 2 . Again using 
Eq. (|J) for standard gravity (without quadratic and non- 
local terms), as a final result one obtains that 



< -(6 2 -3A) 



0, 



(11) 



so that the shear of the homogeneous hypersurfaces 
rapidly approaches zero. One can generalize this result 
to Bianchi type-IX models provided that A is sufficiently 
large compared to spatial-curvature terms || . 

If one is to extend Wald's result to the brane world sce- 
nario, one must consider two new (additional) constraints 
that, taken together, can be seen as sufficient conditions 
for its validity. These new conditions come from the pres- 
ence of the local quadratic term and the nonlocal term 
which, in the early universe, when the isotropization is 
supposed to take place, may play an important role. In- 
deed, when 



S^u^u" < 



(12) 



and 



S^u" > and IS^ - -g^S j u^u u > 0, (13) 

the r.h.s. of Eq. (||) is positive and the r.h.s. of Eq. (0) is 

negative, and Wald's theorem remains valid. 

Let us now separately discuss these two conditions and 
the consequence of their possible violations. The first suf- 
ficient condition imposes a constraint on the projection 
of the bulk Weyl tensor. Using the symmetry properties 
of £ M „, we can decompose it irreducibly as 
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V^u + 2Q( M u I/ ) 



(14) 



where U, <2 M , and are the effective nonlocal energy 
density, flux, and anisotropic stress arising from the free 
gravitational field in the bulk, respectively. The condi- 
tion therefore becomes that 
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U > 0. 



(15) 



Note that this imposes restrictions on the dynamics of 
the metric off the brane. As pointed out a number of 
times, there are no local constraints which can fix the 
sign of U, and this remains a crucial problem for gen- 
eral issues |l4|,[l5| • However it was recently shown from a 
global five-dimensional point of view that for all possible 
homogeneous and isotropic cosmological solutions on the 
brane, the bulk spacetime is Schwarzchild-AdS with the 
mass parameter of the black hole proportional to Ua . In 
this setting one has some constraints on hi {e.g. hi is pos- 
itive for closed S-AdS spacetime p5 16|). For anisotropic 
models such constraints do not exist. 

We therefore choose to consider the possibility that 
U assumes whatever value and sign and we explore all 
possible brane solutions as a function of its value. As 
a first approach we neglect the effect of matter on the 
brane: as we will show, it is straightforward to generalize 
and include a cosmological fluid on the brane, but this 
would obscure the simple analysis we undertake here. We 
then consider coupled equations for <7 A1 „cr A " / , 0, and IA 
in order to study the importance of the bulk term. If 
we restrict ourselves to Bianchi class A models we can 
perform a relatively general analysis without concerning 
ourselves with off diagonal terms in the Ricci tensor. For 
this class of models we thus have Q M — 0, which simplifies 
the set of the equations to solve. 

From the nonlocal energy conservation we have 



U + -QU + a^V^u = 0, 



(16) 



and from the Gauss-Codazzi equation we have an evolu- 
tion equation for the shear, 

d-ftu + Qo-^ + ^R^ - = Q TV; (17) 

where ^'R^u is the Ricci tensor of the homogeneous hy- 
persurfaces. We first nevertheless restrict ourselves to 
models with ^'R^ v — 0, e.g., Bianchi type-I models. 

The system of equation obtained by combining Eqs. 
(§)> (0). ©! and © is not closed. This is a par- 
ticular example of a more general result: due to the 
lack of an evolution equation for the nonlocal anisotropic 
stress Vfiv, the evolution of physical quantity on the 
brane cannot be predicted without ma king further as- 
sumptions on the bulk configuration fl^Jllfl . We there- 
fore choose to set a^V^ = 0. As already noted in ||, 
this condition implies an evolution equation for V^v, i.e., 
V^o^ + {k/k^V^V^ = 0, which can be derived from 
Eq. (^) and is consistent on the brane, allowing us to 
close the system. We do not know whether this corre- 
sponds to a plausible physical form of the bulk metric, 
but this assumption is completely consistent from what 
concerns the brane point of view. Other assumptions un- 
der which one can still solve the system are nevertheless 



possible and we will mention some of them at the end of 
the paper. 

Solving for the nonlocal energy density IA and the ex- 
pansion rate O, we can perform a stability analysis of 
the dynamics by defining a new set of dimcnsionless vari- 
ables: 



T = t 



A 



1/2 



, x = 



' fJ,V 







Y = U = 

(3A)V2' 



2A 



(18) 



This yields a second order autonomous system of nonlin- 
ear differential equations (with ' = d/dr), 



X' + 6XY = 0, 

Y' + Y 2 + U - 1 + 2X = 0, 

U' + AUY = 0, 

plus a constraint equation, 

Y 2 -X -U = l, 



(19) 



(20) 



which can be solved to find the asymptotic evolution of 
the W-term, the shear, and the expansion. Note that, 
if one considers the constraint equation with [7 = 0, 
Y 2 — X — 1 = 0, one recovers the evolution equations 
of standard gravity for the shear and expansion rate. 

This system of equations is quite simple to analyze, 
and reveals the novel behavior of the asymptotic evo- 
lution of the brane. One can identify three critical 
points: (X, Y,U) = (0,1,0), which is an attractor and 
is the same fixed asymptotically stable point we find 
in standard four dimensional cosmology (0 — > (3A) 1 / 2 , 
tr^cTpy -> 0); (X,Y,U) = (0,-1,0), which is a re- 
seller, an unstable point in the contracting phase; and 
(X, Y, U) — (2, 0, —3), which is a saddle point, and there- 
fore unstable. The latter is the most interesting point of 
the system since it corresponds to = and around this 
point the brane "decides" whether to enter the de Sitter 
phase or to collapse. 

The orbits for this system of equations are given by 



X = c\u\ 3/2 , 

u + c\u\ 3/2 , 



Y 2 = 1 



(21) 



where C is an integration constant which depends on the 
initial conditions. The direction field is shown in Fig. 
1. The thick red line represents the condition X = 0. 
Our physical phase space resides where the shear is non- 
negative, i.e., Y 2 > l + U. The black line passing through 
the saddle point represents, instead, the separatrix be- 
tween the region of asymptotic stability (above) and the 
one corresponding to instability, and is given by 



Y 



-(l + [/ + C |[/| 3/2 ) 1/2 , U<-3 
(1 + C/ + Colt/I 3 / 2 ) 1 / 2 , U>-3, 



(22) 
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FIG. 1. Direction field in (U, Y) phase space. The thick 
(red) line represents the condition X = 0. The thin (black) 
line separates the region of asymptotic stability from the re- 
gion of instability of the brane. For values of U < —3 the 
brane may collapse even if it isinitially expanding. 

where Co = 2/3 3 / 2 . The dynamics is therefore such 
that, for any initial condition above this line, the sys- 
tem asymptotically converges to (X = 0, Y = 1, U = 0), 



which corresponds to a^a^ = 0, 8 = V3A, and U = 0; 
however, for initial conditions below this separatrix, the 
brane universe collapses and the anisotropy grows to in- 
finity. 

Isotropy is reached for any initial condition correspond- 
ing to U > 0. We therefore recover the sufficient condi- 
tion for the validity of Wald's theorem mentioned already 
above. When U < the brane may not isotropize and 
can instead collapse. The conventional four dimensional 
cosmology resides on the subspace (the Y axis) given by 
U = 0. As a side result, we can see that when the brane 
enters the de Sitter phase, the IA term approaches zero 
asymptotically. Furthermore, there is a region of initial 
conditions corresponding to small hi where it is possible 
to start with the brane contracting and still reach the sta- 
ble point. This reflects the new properties of the modified 
gravity and the possibility of avoiding the cosmological 
singularity in the presence of the nonlocal term U, which 
was already noted in Q . 

What we have found is not really surprising: the 
projected bulk Weyl tensor plays a crucial role in the 
isotropization of the brane. This means that there are 
additional constraints one must impose in the bulk if one 
is to find a complete generalisation of Wald's theorem to 
the brane; in some sense this defeats the purpose of such a 
no-hair theorem. However, we have been able to present 
a local analysis in which we have closed the system by 



simplifying the effect of the bulk and still obtained non- 
trivial dynamics. Note that we can obtain trajectories 
with U < 0; however, this is not the crucial diagnostic 
for whether the brane isotropizes or not. Of course it is 
of paramount importance to see what global constraints 
the bulk imposes on U in the same way as has been done 
for homogeneous and isotropic models. 

Let us go back to the second of the sufficient conditions, 
Eq. ( |l3| ) . We can covariantly decompose the brane energy 
momentum tensor T^ u into its irreducible decomposition, 



(23) 



where p and p are the energy density and isotropic pres- 
sure, is the energy flux, and 7r M „ is the anisotropic 
stress of the matter on the brane. 

The restrictions for the quadratic term therefore be- 
comes 

2p 2 - 371-^71-^ > and 2p 2 + 3pp - 3<^ > 0. (24) 

These conditions are trivially satisfied for the case of a 
perfect fluid (which satisfies the energy conditions) and 
for a free homogeneous scalar field with a velocity flow 
coincident with u^. However they may not be satisfied 
in general. 

An unambiguous example of violation of one of the 
above conditions is given by a cosmological pure magnetic 
field -B M , such that B^u^ = (for an extensive list of 
anisotropic fluids and sources see [0); in this case we 
have that 



2pj 



z*W u = ~W< o, 



where 



and 



2p 2 B + 3p BPB - = -B 4 > 0, 



Pb = 'Zpb = B 2 /2 



B^By 



(1/3)V£ 2 



(25) 



(26) 



For simplicity we restrict ourselves to the case of 
Bianchi type-I cosmology, which can naturally accommo- 
date a spatially homogeneous ordered magnetic field, and 
we neglect any other form of matter. The presence of the 
magnetic field introduces a preferred direction and we can 
consider the possibility that the background anisotropy 
is due to this presence. This corresponds to the typical 
situation in which the magnetic field is a shear eigenvec- 
tor, 



UauB^ = cr z B„, 



(27) 



which does not represent the most general case; however 
it considerably simplifies the equations . Wc consider 
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the back-reaction from the bulk but, to be able to close 
the system, we assume a^ u P^ u — and P^B^ = 0. In 
this case, together with the Raychauduri equation, the 
full set of equations to solve becomes 

PB+e(pB+PB)+<J^n^ = (28) 

from the local energy conservation, 

4 k 4 

U + -QU = —[6n^+6(p B +p B )a^Tr^ 

+ 2Wir ia> - 2a^TT afl n a v ] (29) 
from the nonlocal conservation equation, and 

k 4 

o><7^ + Qa^a^ = K 2 -k^ v - ~^[{pb + Spb)^^ 

- tt^tt/]^ (30) 

from the Gauss-Codazzi equation. 
Defining the additional new variables 

/ o \ 1/2 

W = a z (jj and Z= ^ (31) 

we obtain an autonomous nonlinear system, 

X' + 6YX + A X W[Z - (4/3)Z 2 ] = 0, 
Y' + Y 2 - 1 + 2X + 2 X [Z + (3/2)Z 2 ] + U = 0, 
Z' + AZY - 2WZ = 0, (32) 
U' + AUY - A x Z 2 [{b/Z)W - 6F] = 0, 
W + 3YW + 8 X [Z - [A/3)Z 2 ] = 0, 

together with the constraint equation 

Y 2 -2 X [Z - {3/2)Z 2 ]-U - X = 1. (33) 

Here X = k 2 X/2A. 

This system has at least three critical points which 
correspond to those obtained analyzing the case without 
magnetic field: (X = 0, Y = 1, Z = 0, U = 0, W = 0) 
is the de Sitter attractor, (X = 0,Y = -1, Z = 0,U = 
0, W = 0) is the repeller, and {X = 2, Y = 0, Z = 0, U = 
—3) is a saddle point for any value of W. They are all at 
Z = and, indeed, p B = is an invariant submanifold 
of the system, i.e., if an orbit starts with ps — it will 
maintain this condition. When the magnetic field is much 
weaker than the brane tension Z 1, we can neglect all 
the Z 2 contributions. In this case we still find a region of 
initial conditions corresponding to collapse of the brane, 
but only for U < —3; this means that the negativity of 
the U-teim is responsible for the collapse and not the 
magnetic field. 

When neglecting the Z 2 contribution we can eventu- 
ally set U — and recover standard gravity. Note that 
this is not possible when the quadratic term is included. 
Indeed, by analyzing the full system of differential equa- 
tions one finds that U = is not a consistent solution. 



By setting U and its derivative to zero we have a conser- 
vation equation involving quadratic terms of the energy 
momentum tensor and implying W — (18/5)1", which 
is not consistent with the evolution equations for Y and 
W. We therefore find an interesting result: in this setup 
we cannot set the nonlocal energy density to zero unless 
we neglect quadratic terms. It is conceivable that this 
is a consequence of discarding the effect of the nonlo- 
cal anisotropic stress, and we are aware that without a 
full analysis of the bulk to completely solve the evolution 
equations this must remain an open problem in this work. 
This example, nevertheless, underline the importance of 
the back-reaction of nonlocal bulk effects in the presence 
of an anisotropic stress on the brane. (As a perturba- 
tion, the effect of anisotropic stresses was considered and 
discussed in p9[.) 

When X = the evolutions of X, Y, and U are decou- 
pled from Z and W, and the phase space is given by the 
product (Z,W) x (X,Y,U)z=o- The existence of other 
critical points depends on the value of the parameter X . 
For 

< X < ^, (34) 

there are three more saddle points, one at (X = 2 — 
(27/8) X , Y = 0, Z = 3/4, U = -3 + (57/16)*, W = 0), 
and other two at Y ^ 0. These saddle points degenerate 
to one when X — 16/27 at the bifurcation point (X — 
0,y = 0,Z = 3/4, U = -3 + (57/27), W = 0). For 
X > 16/27 there is no other saddle point. 

The novelty in the behavior of this system is the fact 
that there are regions of phase space where the brane is 
unstable even when the U-teim is positive. These regions 
can be found by numerically integrating system (|32] ) and 
they typically appear for large value of the magnetic field 
Pb ^ A. Therefore the introduction of the anisotropic 
stress present in the quadratic term of the projected Ein- 
stein's equations makes the brane unstable even when the 
U > condition is satisfied. We conclude that the viola- 
tion of the second sufficient condition for the validity of 
Wald's theorem is enough to cause an expanding brane 
to collapse, provided that the source of violation is suffi- 
ciently strong. 

In order to visualize the five-dimensional phase space 
and have an idea of the orbits of the system we use the 
fact that W 2 = AX is another invariant submanifold. It 
corresponds to the axial symmetry of the brane universe 
around the z direction and is realized when a^ v oc tt^. 
This condition allows us to get rid of the evolution equa- 
tion for X and use the constraint equation to eliminate 
U. We are thus left with a closed system of three dif- 
ferential equations in Z, W, and Y. In Fig. 2 we plot 
some typical trajectories of this system for X = 1/2. As 
one can see, for large value of the energy density of the 
magnetic field, Z > 1, the brane can collapse. 

Let us discuss some possible generalizations of our re- 
sults. In the analysis undertaken so far, we have assumed 
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FIG. 2. Phase space corresponding to the invariant sub- 
manifold W 2 = AX. For large values of Z, typically Z > 1, 
the brane can collapse, even if U is positive. Note that, since 
W 2 = 4X, this phase space cannot contain the new saddle 
point at Z — 3/4. 

that the curvature of the homogeneous spatial hypersur- 
faces is zero. However one can still generalise some of the 
results obtained for the case ps = 0, by including the cur- 
vature term in the system of equations (|2^) . Defining the 
dimcnsionless variable 



the relevant equations become 

Y' + Y 2 + U -1 + 2X = 0, (36) 
U' + AUY = 0, 

with the constraint 

Y 2 -X- U- 1 = -R/2 > 0, (37) 

where we have used that ' 3 'i? < 0. This, as already 
stated, is true for all Bianchi models except type-IX. 

Let us define with Yr—q (U) the "old" orbits of the 
system obtained by setting the curvature term to zero, 
and given by Eq. By carefully integrating the evo- 

lution equations for Y and U and using the inequality 
(p7|) , one obtains the following behavior of the "new" or- 
bits Y(U): if at some point of the phase space of Fig. 1 
Y(U(0)) — Yr = o(U(0)) (a new orbit crosses an old one 
at r = 0), then 

Y(U(t)) > Y r=0 (U(t)) } r > 0, 
Y(U(t)) < Y r= o(U(t)), t < 0. 

By again considering Fig. 1 and applying these inequali- 
ties, we can easily conclude that the region of asymptotic 
stability is enlarged by the introduction of the curvature 



term. Note that the orbits never cross the Y axis since 
there U' = 0. Indeed we still find the attractor and the 
repeller at the usual positions (X, Y, U) = (0, ±1, 0), cor- 
responding to R = 0. However, it is impossible to show 
the existence of a saddle point at Y = unless one can 
show that R is stationary at this point. This is true for 
models in which the surfaces of transitivity are isotropic 
(see pi), 

(3 H- = I (3) Rh^, (38) 

but is not generally the case. Looking at Fig. 1, in these 
models the saddle point is displaced towards the left, at 
U = -3 + R. 

According to these inequalities it is still possible that 
the universe recollapses. Other types of non de Sitter 
behavior are of course nevertheless possible, in particu- 
lar, an asymptotic approach to an Einstein static brane. 
However, even with the introduction of a curvature term, 
it seems reasonable to conjecture that the only types of 
stable asymptotic behaviors are an asymptotic approach 
to de Sitter spacetime and recollapse, as in the case of 
^Rfj.u = 0. It is well known that the Kasner (vacuum 
Bianchi type-I model) collapsing solution is an attractor 
of Bianchi class A models with ordinary matter content. 
It would of course be important to understand if this is 
still the case in the presence of a nonlocal bulk effect. We 
reserve this study for future work. 

The effect of matter and other possible cosmological 
fluids has been discarded in our analysis, but it can be 
included quite straightforwardly by adding a perfect fluid 
moving orthogonally to the homogeneous hypersurfaces 
with state equation p = ujp to our system of equations. 
The presence of a new variable p introduces an extra 
dimension in the phase space and a new parameter u>. 

When the magnetic field is not taken into account we 
can integrate the orbits of the system using the dimen- 
sionless variable Zp = p/X. This yields 

x = c\u\ 3/2 , 

Z F = D\lJf" +1 V\ 
Y 2 = 1 + U + C\U\ 3/2 

+ x {w\u? { " +1)/i + £ 2 |t/| 3(w+1)/4 ), 

where D is a new integration constant. The attractor 
and the repeller are at the usual positions (X,Y,U) — 
(0, ±1, 0) and Zp — 0. However, now the saddle point at 
Y = becomes a line which depends on Zp: 

U = -3(1 - x(w - l)Z F - X uZ%). (39) 

Note, however, that for lo > —1/3 the allowed values of 
Zp of this line are limited by the constraint X > 0. The 
separatrix between the region of asymptotic stability and 
instability is the two dimensional surface which forms the 
saddle line at the intersection with Y = 0. The region of 
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instability always lies below this separatrix. The general 
characteristics of the phase space of the system, as illus- 
trated in Fig. 1, remain unchanged: the only two stable 
asymptotic behaviors are an approach to the de Sitter 
spacetime and recollapse; recollapse can only occur when 
U is negative. On the U = plane, we recover the results 
of 0. 

In order to close the system of equations discussed 
above, we have always chosen to neglect the nonlocal 
anisotropic stress V^. As already discussed, this always 
corresponds to a consistent solution in the bulk, pro- 
vided that the equations are consistent on the brane, al- 
though it may not correspond to some plausible physical 
bulk configuration. For future investigations we mention 
three other possibilities which allow one to consistently 
close the system on the brane. If you have some sort of 
anisotropic stress on the brane a very simple possibility 
is of course to assume that U = U = 0. In this case 
one obtainss a consistent closed system provided that a 
nonzero V^ v is considered. 

A somewhat more nontrivial assumption is the fol- 
lowing. Since the nonlocal anisotropic stress Vp, v is the 
spatial component of the projected five-dimensional bulk 
Weyl tensor on the brane, it is physically plausible to 
have a proportionality between this term and the electric 
part of the brane Weyl tensor, 

a(^\ V l tv=E^ = C mv pu a u , (40) 

where C^ av /s is the Weyl tensor on the brane. With the 
help of the spatial traceless component of Einstein's equa- 
tions [see pl|| , Eq. (A3)] one can find E^ u in terms of the 
shear and the Ricci tensor of the homogeneous hypersur- 
faces and have an expression for V^u which depends on 
a constant parameter a. 

Another possibility to explore is to assume that V^v 
evolves according to a "free" evolution equation, 

7V + 067V = °> ( 41 ) 

where j3 is a constant parameter. In both cases V^v feeds 
into the propagation equation for the shear and the 14- 
tcrm, possibly giving rise to interesting dynamics, which 
can be studied by varying a and (3. 

We briefly summarize the present work. We have de- 
rived a set of sufficient conditions to be satisfied by the 
local brane matter and the nonlocal energy density, so 
that a homogeneous and anisotropic brane with tension 
A, satisfying modified Einstein's equations with a positive 
effective cosmological constant, asymptotically evolves to 
de Sitter spacetime. We have discussed the violation 
of these conditions, showing that, in the presence of a 
negative nonlocal energy density hi < and/or strong 
anisotropic stress \it^ v \ > A, a Bianchi type-I brane, even 
if initially expanding, is unstable and may collapse. This 
is our main result. We have also given an explicit example 



of source of anisotropic stress by considering the case of 
a cosmological magnetic field on the brane, and we have 
shown that the back-reaction of nonlocal effects cannot 
be neglected. We have discussed possible generalisations. 

Finally, let us reemphasise a key aspect of our anal- 
ysis. We have studied the asymptotic behavior of the 
brane from a local point of view. On the one hand, this 
approach is incomplete and this is evident in what we 
have found: the bulk back-reaction plays a crucial role in 
the evolution of the brane. This back-reaction has been 
discarded in previous local analyses where isotropic flu- 
ids were considered [this is in fact consistent with our 
analysis, merely corresponding to a choice of integration 
constant when solving Eqs. (^0|)]. Given that we do not 
have an evolution equation for the nonlocal anisotropic 
stress Vfivt w e have only considered a subspace of pos- 
sible solutions such that we could close the system. We 
should point out, however that the complementary ap- 
proach, i.e., constructing complete bulk solutions and as- 
sessing their effect on the brane, may also be restrictive 
although the two approaches are clearly necessary for a 
complete understanding of the asymptotic stability of the 
brane. 
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